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Abstract 

I> | We consider the renormalization of theories with many scalar fields. We discuss at the one- 

loop level some simple, non-gauge models with an arbitrary number of scalars and fermions 
both in mass-shell and MS schemes. In MS scheme we give a detailed qualitative analysis of 
iy-j ■ the RG flow of dimensionless couplings in flavor space. 
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1 1 Introduction 



Q_i' In this paper we consider the renormalization of theories with many scalar fields which mix non- 



trivially under interactions. Specifically, we discuss at the one-loop level two simple, non-gauge 



models involving scalar fields and fermions. We understand mixing of fields in a broad sense, as 
the existence of flavor off-diagonal propagators that do not vanish. (The stronger condition that 
they should not vanish on the mass shell is necessary for mixing of mass eigenstates, but it will not 
be required in what follows.) Mixing is trivial if a flavor basis exists such that there is no mixing in 
that basis. We refer only to orthonormal flavor bases in which the kinetic terms are diagonal and 
normalized, and to unitary or orthogonal transformations among those bases. Therefore, mixing 
is trivial if the relevant two-point function can be diagonalized by a (coordinate and momentum 
independent) unitary transformation, i.e., if it is normal as a matrix in flavor space. In the case of 
fermion fields renormalizable, Lorentz-invariant interaction terms can only be bilinear in the fields. 
In the scalar sector a larger variety of interactions is allowed and, as a consequence, mixing can be 
due to interactions that are quadratic, cubic or quartic in the scalar fields. 

In many theories, like the Standard Model, there exists an "interaction" flavor basis which 
diagonalizes the bilinear interaction terms. More generally, we define an interaction basis as one 
in which there would be no mixing if mass matrices were proportional to the identity. Such bases 
need not exist in general, however, as would happen in a model with two different fermion currents 
which are not simultaneously diagonalizable. This phenomenon is even more obvious in theories 
with mixing of scalar fields. Due to the presence of the different types of interactions mentioned 
above it is often impossible to find a flavor basis in which there is no mixing, even if mass terms 
are trivial. Our treatment is independent of the existence of an interaction basis. 

In the next section we study a model with N s scalar and Nf fermion fields. The scalar fields 
mix through a bilinear term involving a mixing matrix (and an additional field taken also to be a 
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scalar), quartic self-interactions and Yukawa interactions. The model does not possess continuous 
symmetries and is intended to minimize the effect on mixing of the quartic terms in the potential 
as much as possible. We discuss the renormalization of this model in detail at the one-loop level, 
for generic N s and Nf, both in on-shell (OS) and MS schemes. In MS we pay particular attention 
to the renormalization group (RG) flow of dimesionless parameters. 

The treatment of counterterms in section 2 follows closely our previous paper [1] on fermion 
mixing, with which we try to keep the overlap at a minimum. Besides the different nature of the 
interaction Lagrangian for scalar and fermion fields, another important difference with [1] is that in 
this paper we consider the case of several coupling matrices, Yukawa couplings being row matrices, 
whereas [1] is restricted to models with only one. The presence of several coupling matrices is not 
only algebraically more involved but also qualitatively more interesting as reflected, for instance, 
in the RG equations. 

In section 3 we consider a model with spontaneously broken SU(N) symmetry which is, in fact, 
the matter sector of an SU(N) Georgi-Glashow model coupled to fermions. This is quite different 
from the model in the previous section, since mixing is controlled by the continuous symmetry 
and its patterns of spontaneous and explicit breaking, and it occurs mostly due to the quartic self- 
interactions. We restrict ourselves to MS in that section and to one particular type of spontaneous 
symmetry breaking (SSB). We discuss in some detail the simpler case in which there is no explicit 
SU(N) breaking, and consequently no coupling matrix, and comment briefly on the case with 
the symmetry explicitly broken by mass terms, in which an orthogonal coupling matrix enters the 
Lagrangian. In section 4 we give our final remarks. Complementary material is gathered in several 
appendices. 



2 A model of scalar mixing with N s scalars and Nf fermions 

We consider a schematic model involving N s real scalar fields (p = . . . , 4> Ns ) and Nf fermions 
ip = (if) 1 , . . . , i/j n ) interacting through a Yukawa coupling. The Lagrangian includes an additional 
field x n °t coupled to fermions, chosen to be a scalar for simplicity, C = £ frcc + £ int with 

£ frco = ~x (□ + m 2 x ) X ~ \<\>- (□ + M 2 ) 4> + 

£ A 1 Nf (1) 

Ant = ~y X 3 ~ ^jjX 4 - Vi(<f>) - -<j>-H<p X -J2a k -(f>^ k ?Pk ■ 

fc=i 

The <fi mass matrix M 2 is real symmetric, positive-definite and regular (i.e., none of its eigenvalues 
is degenerate). The mass-degenerate case is discussed separately below. The <p~x interaction term 
is bilinear in <fi, with a real symmetric coupling matrix H. V^(cf)) contains only terms quartic 
in <f), whose explicit form is considered in detail below. Yukawa interactions are specified by Nf 
real coupling vectors a k = (a kl , . . . , a kNs ), k = l,...,Nf. The field <fi takes values in real N s 
dimensional Euclidean space, which we refer to as "(scalar) flavor space." We do not use the 
summation convention for the flavor indices of scalar and fermion fields. 

C is invariant under an exact U{l) N f symmetry of the fermion sector, related to the conservation 
of the number of each fermion species ip k , which therefore do not mix among each other. The 
fermion fields tj) are chosen to be massless. C is then also invariant under the exact Z2 symmetry 
4> — > —0, V ~^ TbV'j which forbids vertices involving an odd number of 4> fields, couplings of x to 
fermions, and x~4> mixing. An interaction term of the form (fi-Kcfrx 2 with a coupling matrix K is 
allowed by the Z2 symmetry, but we do not include it in C for simplicity. 
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Let us call S Y the subspace of flavor space spanned by {a k } k Li aR d d Y its dimension, so d Y < Nf 
and d Y < N s . If d Y < N s , the Yukawa interaction is invariant under the discrete symmetry Z2y 
that obtains by restricting the transformations of the exact Z2 symmetry of C described above to 
the components of <j> lying on S Y . In other words, under %2y only the components of 4> coupled to 
fermions change sign. If d Y < N s — 2, the Yukawa interaction is obviously invariant under the group 
SO(N s — d Y ) acting in the orthogonal complement of S Y . We assume that V^(cj)) is chosen so that 
it possesses the same symmetries as the Yukawa interactions. Thus, both Z2y and SO(N s — d Y ) 
are approximate symmetries of £, softly broken by x~4> interactions (unless S Y is contained in an 
eigenspace of H) and by the 4> mass term (unless S Y is contained in an eigenspace of M 2 ). 

Depending on the geometric configuration of {a k } there can be other, generally approximate 
symmetries of C If S Y = S Yl © S Y2 with S Yl generated by, say, {a^, 1 < j < rai} and orthogonal 
to 5y 2 , generated by {a j ,n 1 < j < Nf}, we have the approximate symmetry groups a^-cf) — > 

—ay4>, ipj — > 75^, if j < ni, with all other fields unmodified, and similarly 7Li Y2 : a k 4> ~^ — « fc -0, 
tp k — ► J 5 ip k , if «i < fc < Afy, all other fields unaffected. Here we denoted = Oi/||ai||. These 
symmetries are softly broken unless S Yl and (Sy 2 are included in eigenspaces of H and M 2 . The 
converse is also true, if Z,2 Yl or Z2y 2 are exact symmetries of C, then S Yl and 5y 2 must be contained 
in eigenspaces of H and M 2 . 

Similar approximate discrete symmetries appear if the {a k } lie on several different orthogonal 
subspaces, and if within each subspace the coupling vectors are collinear. One extreme situation 
is when all of the coupling vectors {a k } are pairwise orthogonal, with Nf < N s . In that case we 

have a symmetry, in general also softly broken by M 2 and H. The other extreme is when 
all coupling vectors are collinear, so that only one component of <f) is coupled to fermions. In 
the special case in which all of the a k are eigenvectors of H there is an interaction basis (if V4 is 
appropriately chosen), and if they are also eigenvectors of M 2 mixing becomes trivial. Also, if out 
of the Nf coupling vectors a k , n (with n < Nf) are equal, and orthogonal to the remaining ones, 
the Lagrangian is invariant under the corresponding SU(n) symmetry of the fermion sector. The 
stability under renormalization of those special geometric configurations of {a k } is of particular 
interest. As expected, all of them are critical points of the renormalization group (henceforth RG) 
flow. We consider these issues in detail below, in MS scheme. 

Despite its simplicity, in general the Lagrangian (1) does not have an interaction basis in which 
C int is diagonal in <f>, except in the special cases mentioned above. This feature is interesting to 
point out because some approaches to mixing renormalization rely on the existence of a flavor basis 
diagonalizing the interactions that is orthogonally, or unitarily, related to the mass basis. In fact, 
we could further generalize our model to have a Yukawa interaction of the form ^2^=1 ^^-a^fpa, 
where {A a } are N s Hermitian Nf x Nf matrices. C in (1) corresponds to the case where all A a 
are simultaneously diagonal. In the general case where not all those matrices commute there is no 
diagonal "weak" basis for the fermion fields either, even though they are massless. 

2.1 Quartic terms in the potential 

The simplest form for the quartic interactions V^cf)) is isotropic in flavor space. That form is not 
enough to make £ renormalizable, so we adopt the following parametrization for V4, 

V 4 (<f>) = fj(^) 2 + ^ Yj 9ij(ai-<t>){a r <t>)(<t>-(t>) + - Yj 9ijki(^-4>){a r 4>){a k -4>){ai-4>). (2) 

' i<j=l ' i<j<k<l=l 
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We assume that g, gij and gijki are such that V^(0) > for all (p. If N s > d Y + 1, V^(0) is invariant 
under the SO(N s — d Y ) transformations leaving all components a,k-(j> invariant. The coupling 
constants g^ and gijki are chosen so that V^{4>) is invariant under all of the discrete symmetries of 
the Yukawa interactions. Thus, in the example given above of two orthogonal subspaces S Yl and 
S Y2 , we must have g^ = = g^i if an odd number of indices take values < n\. To give another 
example, if the are pairwise orthogonal then only the couplings g, gu and gujj with i < j may 
not vanish. 

2.2 Quantum Lagrangian 

Our perturbation expansion parameters are A, a k = ||a fe ||, g, g^, g^^i and the dimensionful £ and 
H. We introduce % wave-function and mass renormalization constants 5Z X and 5m^ in the usual 
manner. Let us assume for the moment that an interaction term C K = —(fi-Kcfrx 2 is added to £ int , 
with K positive definite. It is easily seen that K must be multiplicatively renormalizable. In MS 
scheme the counterterms to K must be of the form 




with e = 4 — d. The second term in (3) collects divergent e-pole contributions to 5K that do not 
vanish if K = 0. Since 8Z^ K)I originates from diagrams with two <fi- and two \~ external lines and 
no factor of K, it must vanish when H = 0. But 5Z {K)I cannot depend on H [2], therefore it 
vanishes. Thus, if we set K = as done in (1), all Green's functions of <j> a (t>bX 2 are finite once the 
other couplings and the propagators have been renormalized (with K=0). 

There is no big loss of generality in setting K = 0. If we had a term C K in C we would set 
0jv s +i = X an d redefine V4 — > V4 — C K . We would still have the exact Z2 symmetry forbidding 
additional cubic interactions and x~<fr mixing. We would require, however, that C K have the same 
discrete symmetries as the Yukawa term. 

Having set K = 0, a similar argument goes through for A. We can then just set A = 
in C without spoiling its renormalizability. If we also set £ = 0, Green's functions of x 3 receive 
contributions only from triangle-type diagrams with three 4>-Hcj>x external vertices. Such diagrams 
have negative superficial degree of divergence and are therefore finite in the renormalized theory 
with £ = 0. We arrive in this way at a minimal \ sector, which still serves its purpose of inducing 
bilinear mixing of <fi. 

The fermion fields do not mix, so their wave-function renormalization matrix is diagonal, ifj kQ = 

1 /2 

Z^ k ip k = (1 + l/25Z 4 , k )t/) k . The renormalized scalar field (f) can be related to the bare one as 
4> = A(f), with A a real N s x N s matrix. It will prove convenient to introduce the real polar 
decomposition of A explicitly [1], 

4> = QZ*4>, Q = e SQ , Z^ = l + ±SZ, (4) 

where Q is orthogonal, Z symmetric and positive definite, SQ = —SQ l and 5Z = SZ l . In MS 
scheme we set SQ = by definition and, owing to the softly broken 2,2y symmetry, SZ must have 
the form, 

8Z = ^5Z 4>k a k ®a k + ^ 8Z^ k (a,- <g> a k + a k <g> a,) + SZ', with SZ'a t = 0, i < Nf. (5) 

fc=l j<k=l 
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Furthermore, if N s > d Y + 1, due to the softly broken SO(N s — d Y ) symmetry, SZ' must be 
diagonal. In the superrenormalizable case a k = = V4 we must also have SZ = in MS, since 
there is no divergent field renormalization, but even in that case SQ is divergent in OS scheme due 
to infinite mass renormalization. 

We write the relation between bare and renormalized <fi mass matrices as [1], 

M 2 a = Q m t (M 2 + SM 2 )Q m , Q m = e- SQ ™, [SM 2 ,M 2 ]=0, (6) 

with Q m orthogonal, SQm 1 = —SQ m , and (SM 2 ) 1 = 8M 2 . At one-loop level we have M 2 = 
M 2 + SM 2 + [SQ m , SM 2 ], with SQ m determined up to addition of an antisymmetric matrix 
commuting with M 2 . Q m is then defined up to multiplication by an element of the subgroup 
of SO(N s ) that leaves M 2 invariant. We write the counterterms to H similarly [1], taking into 
account that it must be multiplicatively renormalizable, 

H^^V'ZhHV, V = e~ sv , [Z H ,H}=0, (7) 

where V is orthogonal, SV = —SV*, Zh = 1 + SZh = Zh '■ V depends not only on H and the 
coupling constants g, gij, gijki, but also on the unit vectors a k , through V4. For this reason, we do 
not expect to be able to set SV = 0, not even in MS scheme, unlike the case in [1]. At one loop 
we have H = [SV, H] + SZhH, with SV determined up to addition of an antisymmetric matrix 
commuting with H. 

The counterterm parametrizations (6), (7) are a restriction to real matrices of the one given 
in [1] for complex normal matrices. As such, it is applicable in general to real normal matrices 
(i.e., matrices A such that [A, A ] = 0.) For such matrices a multiplicative transformation by 
congruence A = B t AB is not possible if we want both A and its transformed A to be normal, 
except in the particular case when both are symmetric. On the other hand, if A and A are normal 
but not symmetric, it may not be convenient to split them in their symmetric and antisymmetric 
parts and to treat them separately, as is clearly the case for orthogonal coupling matrices. 

Yukawa couplings are multiplicatively renormalizable. Since fermions do not mix, if we set 
some a n = the corresponding fermion field if) n completely decouples. The relation between bare 
and renormalized coupling vectors is of the form aj = (i e / 2 Ajdj, with Aj a real iV s x N s matrix 
depending on {a^} and, in principle, all other couplings in C Any such relation between two 
vectors can be written as the product of a rotation taking aj into a J0 , times a dilatation changing 
aj to cij , 

a kQ =^ 2 Z ak Q ak a k , Z ak >0, Q ak =e~ SQ ^, SQ (lk = -SQ t ak . (8) 

This separate treatment of the renormalization of a k = \ \a k \\ and 0^ = 0^/ a k is necessary for the 
renormalizaton of the potential and turns out to be convenient for the derivation of RG equations. 
At one loop, a ko = ^^(a^ + Sa^) with Sa k = 5Z ak a k — SQ ak a k . The softly broken Wj2 Y symmetry 
of C prevents the subspace S Y from receiving infinite renormalization. Thus, in MS we expect the 
counterterms 8a k to be linear combinations of {a{\. Geometric configurations giving rise to further 
approximate symmetries, as described above, are also preserved in this scheme. In OS scheme the 
symmetry breaking due to x~4> interactions and cj) mass term is apparent in the counterterms, with 
Sa k having finite components orthogonal to S Y . 

We consider, finally, the renormalization of quartic <f> self-couplings. The coupling constants 
g, etc., in V4 receive infinite additive renormalization at one loop from fermion box diagrams. We 
set, 

90 = v e (g + Sg), 9ij = v e (9ij + Sgij), 9ijki = ^(am + Sgijid)- (9) 
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The complete set of counterterms to V 4 is obtained by expanding 



V 4 (0 O , {a ko },g , {9ij }, {9ijki }) = 

V 4 ^QZ 1 / 2 cf>,{Q ak a k }, f i%g + 5g),{fi^g ij +Sg ij )},{fi e (g ijM + Sg^ (10) 

in powers of the coupling constants (see appendix A.l). 

We interpret the classical Lagrangian C in (1) as being written in terms of bare fields, mass and 
coupling constants, and substitute in it the above expressions for bare quantities. The resulting 
Lagrangian, expressed in terms of renormalized quantities, does not depend on the orthogonal ma- 
trices Q, Qm, V and {Q ak } independently, but only on Q m Q, VQ and {Qa k Q}- This ambiguity 
is fixed by imposing appropriate renormalization conditions. In MS and similar schemes we set 
Q = 1, whereas in OS scheme we choose to parametrize the theory in such a way that M 2 {) and 
M 2 are simultaneously diagonal by setting Q m = 1. 

The Lagrangian is then given in terms of renormalized fields and parameters by (1) and (2), 
with A = = £ for simplicity, and with the substitutions H — ► [i e l 2 H, a k — ► [i e a k , g — > //<?, etc., 
plus the counterterm Lagrangian £ ct . Retaining only those counterterms that contribute at the 
one-loop level, C ct has the form, 

C ct = - l -8Z xX u x - l -5m 2 xX 2 - \<t>-8ZU<t> - ^AM 2 

Nt N 



+ Y i SZ k ^ k i^ k - i/// 2 0-Aif0 x - v e/2 J2 Aa k -<W k il> k + SV 4 (11) 



with, 



k=i fc=i 



AM 2 = SM 2 + [SQ + SQ m , M 2 } + ^{SZ, M 2 }, 
AH = [SV + SQ, H} + ± {SZ, H} + SZ H H + X -SZ X H, (12) 
Aa k = (dZ ak + SZ^ - SQ ak + SQ + ±6z\ a k , 



and SV 4 as obtained from (10) in appendix A.l. 



2.3 On-shell scheme 

In OS scheme we set SQ m = 0, and require M 2 to be diagonal to all orders in perturbation theory. 
We assume M 2 to be regular and then, since [SM 2 , M 2 ] = by definition, SM 2 and M 2 must 
both be diagonal. We assume for simplicity that \ an d <p masses are such that \ ls stable, and 
that decays 4> a — > (f>b + X are n °t allowed. The masses of stable particles, fermions and %, are pole 
masses in this scheme. 4> states are unstable. We adopt the OS renormalization scheme for the <fi 
propagator as well (see e.g. [3]). 

The x field 2-point function T2 X (p 2 ) = p 2 — m x — H x (p 2 ) is obtained by imposing the renor- 
malization conditions Ii x {m x ) = = Ii' x (m x ), where the prime stands for d/dp 2 . Denoting £l x {p 2 ) 
the unrenormalized x self-energy at one loop, 

1 1 Ns 

^ p2) = ~ 16^ Tr ^ + 32^ E H ab H ba b (p 2 ,m 2 a ,m 2 ), (13) 

a,b=l 
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we get, 

5Z X = f^.(m 2 ), <5ra 2 = m x^'x^ m 'y) ~ ^x( m x) (14) 

r 2x ( P 2 ) = (i + n' x {m x )) ( P 2 - m x ) - (n x ( P 2 ) - n x (m x )) , (15) 

with tt' x (m x ) < if x is stable, and T2 X (p 2 ) independent of the dimensional regularization scale 
fx. In (13), 60 refers to the finite part of the 2-point loop integral, as defined in appendix B. As 
expected, SZ X is finite. 
The 2-point function 

r 2 ( P 2 ) = p 2 i-m 2 -n( P 2 ) 

in this scheme is obtained with the renormalization conditions, 

mtU aa {m 2 a ) = = ^tW aa (m 2 a ), ^tU ab (m 2 a ) = = mtU ab (m 2 b ) if a ^ b. (16) 

Here, the renormalized self-energy function H(p 2 ) = ^l(p 2 ) — p 2 SZ + AAf 2 is given in terms of 
the unrenormalized one, 

si(p 2 ) = rt 1 \ P 2 ) + nW(p 2 ) + rt :i \ 



n £V) = -^ e (H 2 ) ab + ^J^H an H nb b (p 2 ,rn 2 x ,m 2 n ), 

n=l 

j N f 1 N f 

OP) (p 2 ) = -—^-p 2 a k ®ak + ^P 2 b (p 2 ,0, 0)Y, a k® a k 



(17) 



4vr 2 e^ 1 8vr2 ; 

k=l k=l 



where 

b (p 2 , 0, 0) = log -18^-2 = log (|i ) - i7rQ(p 2 ) - 2, 

the second equality being valid for real values of p 2 . The term fl^ in (17) gathers the contributions 
to f2(p 2 ) from diagrams with one insertion of V4. These are independent of p 2 , real, and symmetric 
with respect to flavor indices (we give their explicit expressions for the case N s = 2 in appendix 
A). Cl(p 2 ) is then symmetric in flavor space. 

Applying the renormalization conditions (16) we get, for the diagonal self-energies, 

IWp 2 ) = n aa (p 2 ) - men aa (m 2 a ) - (p 2 - m 2 )<neO' aa (m 2 ), (18a) 
5Z aa = mcn' aa (m 2 a ), AM 2 a = m 2 a ^n' a M 2 a)-^aa(m 2 a ), (18b) 

and for the off-diagonal ones, 07^6, 

99 99 

IU(p 2 ) = tt ab (p 2 ) - \y\tSl ab (m 2 a ) + \ \KtSl ah (m\\ (18c) 

m£ -ml m z a - m z h 

6Zab = Ktn ab (m 2 a ) - fttn ab (m 2 ) = m 2 V\t?l ah {m 2 a ) - m 2 a ^n ab (m 2 ) 

m 2 a -m 2 b ' ab m 2 a -ml 

One-loop diagrams with a V4 vertex do not contribute to wave-function renormalization, so that 
SZ is independent of ft^\ Equation (18) then gives an explicit expression for SZ, which is of 
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course symmetric, and finite if all a k = 0. The expression for AM 2 depends on the detailed form 
of fi (3) . From (12) with SQ m = we get 

SM 2 + [SQ, M 2 } = AM 2 - - {SZ, M 2 } . 

Substituting the counterterms on the r.h.s. by their values from (18), and taking into account that 
SM 2 is diagonal in this scheme and SQ antisymmetric, we obtain, 

8M 2 a = -Xtn aa (m 2 a ), 5M 2 b = if a ± b, 

xn 1 Ktn ab (m 2 a ) + dltn ah (m 2 b ) (19) 
2 m~ — m b 

Notice that SQ is divergent even if all a k = (i.e., ft (2) = 0) and V 4 = (i.e., ft (3) = 0), due to 
divergent terms in fi^, although in that case the theory is superrenormalizable. 

The expression (19) for SQ is valid only for regular mass matrices. If some of the squared 
masses m 2 are degenerate the derivation given above for SQ does not hold, since in that case 
the renormalization conditions (16) for off-diagonal 2-point functions are not independent. The 
treatment of the renormalization of 2-point functions T2 in the mass-degenerate case is completely 
analogous to that given in [1] in the case of fermion fields. 

The fermion two-point functions have the form T2^ k (p) p / H- l p k (p 2 ), with 



1 1 2 

n^ fc {pf = SZ^ k - J^2^ a l + 327r2 2 Yl {ak)c 2 (m 2 c a (m 2 c ) + (p 2 - m 2 c )b (p 2 , m 2 , 0)). (20) 

c=l 

The renormalization condition T 2l p k (p)u k (p) = if p 2 = is trivially satisfied, whereas the condition 
l/V r 2y> fc (pW(p) = u k [p) \ip 2 = 0, or n^ fc (p 2 = 0) = 0, leads to, 



-1 1 1^2 

SZ *>> = i^ fl * + + ^2 E (a k )c 2 a (m 2 c ), 



16vr 2 e k 64tt 2 k 32vr 2 

c=1 (21) 



c=l 



m 2 c ){bo{p 2 ,m 2 c ,Q) -a (m 2 )), 



where we used 6o(0, m 2 ,0) = ao(m 2 ), and 6' (0, m 2 ,0) = — l/(2m 2 ), as is easy to check with the 
expressions given in appendix B. 

We consider next the <p~x vertex function, which has the form, 

r ab (pi,P2) = -v e/2 H ab - fi e/2 AH ab + r£(pi,P2) + r^\ Pl , P2 ). (22) 

Here, r^ b \pi,p 2 ) is the contribution from the scalar triangle diagram, 

r a b (Pi,P2) = -M e/2 Y H acH cd H db Co{pi + p 2 , p 2 , m 2 , , m 2 d , m 2 c ) , (23) 

c,d=l 

where Cq is the finite scalar triangle integral, as given in appendix B. There are no other triangle 

(2) 

diagrams at one loop, since we set £ = 0. T ab (pi, p 2 ) collects the contributions from diagrams with 
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one V4 vertex, which are divergent and depend only on p\+p2- Its explicit form in the case N s = 2 

■ ii ) (pi,P2) = r£ ) 



is given in appendix A. It is apparent from (23) that (pi,P2) = ri„ (j>2,Pi), and from appendix 



m 



x: 



(2) (2) 

A that ry (pi,p 2 ) = r^(p 2 ,Pi), therefore, T ab (p 1 ,p 2 ) = T ba (p 2 ,pi). 

If the external momenta in T a b(pi,p2) are on mass shell, pi = m 2 , p\ = m 2 , and (P1+P2 
Defining 

ri? s) = r a6 (pi,p 2 )| on _ sheII , (24) 

we see that T^ S ^ does not depend on momenta, since it can only depend on p\, p\, p\-p2 which 
are fixed by the on-shell conditions, and furthermore ^ = T^ a ^ ■ For concreteness, we impose 
on T a b(pi,p2) the renormalization condition T^ >S ' ) = H ab , whence, 



AH. 



ab 



riiV,P2) + r£V,P2)l . (25) 

J on— shell 



The one-loop counterterm AH is completely fixed by this equation. We can, however, determine 
the counterterms to the coupling matrix as defined in (12) through the equation, 

HSZ H + [SV, H] = AH - [SQ, H] - i {SZ, H} - ^SZ X H (26) 

where all the quantities on the r.h.s. are already known. This equation can always be solved for 
8Zh and SV, as shown in appendix C of [1], by projecting over an appropriate basis of matrices. 
In appendix C we consider the special case in which H is regular. 

We consider, finally, the renormalization of the Yukawa couplings, as given by the 1-PI three- 
point Green's function, T ka (p\,p2). Here, p\ t 2 are the momenta of the fermions, incoming by 
convention, k = l,...,Nf is the fermion flavor, and a = 1,...,N S the scalar flavor. Only one 
fermion flavor index is needed since, due to the U(l) N f symmetry, Green functions of ^aV'fcV'j 
vanish unless j = k. T ka (pi,p2) is a bispinor which at tree level is proportional to the identity 
matrix, T ka (pi,p2) = (a-k)a- Expanding the one-loop T ka (pi,p2) in the usual Dirac matrix basis, it 
is clear that only the scalar form factor can receive divergent contributions, the other form factors 
being finite. We focus then on ^ k S J{pi,P2) = l/4Tr£>rfc a (pi,£>2), assuming for concreteness that 

the physical values of the coupling vectors a k are fixed by these form factors. 

(s) 

At one loop T ka is given by, 

ri?(pi,P2) = S /2 (a k ) a + ^'\Aa k ) a + f' 2 {a k ) a £ (a k ) c 2 ( Pl ,p 2 , m 2 c ), (27) 

c=l 

where is the scalar form factor for the triangle diagram, 

f iS) (pi,P2, ml) = -=1 + (p Q (pl ml 0) + b (p 2 , ml, 0)) - ^±^C ( P i + p 2 ,P2, m 2 c , 0, 0). 



8vr 2 e 32vr 2 

On the mass-shell p\ = = p\, 1p\-p 2 = m 2 . Defining 



(28) 



r£ S) = rf a } ( M ) h , fW° S \m 2 a ,m 2 c ) = f^(Pi,P2,m 2 c ) , (29) 

on— shell on— shell 

and setting 9\tT k ^ = fi e / 2 (a k ) a , from (27) we have, 

(Aa k ) a = -(a k ) a Y,(*k)c 2 Kcf^ os \m 2 a ,m 2 c ). (30) 
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These counterterms are enough to renormalize the scalar form factors at one-loop level. The 
renormalization constants 8Z ak and 8Q ak defined in (8) and (12) can be extracted from the value 
of Aaj.. For that purpose, it is convenient to rewrite /( 5 ')( os ) as, 

mtf^ os \ml,ml) = - JL + /W(m») + fW{ml,m*), (31a) 

= J^Mml), f™(ml,m 2 c ) = ^ WeC„(pi + p 2 ,p 2 , m 2 c , 0, 0)| on _ shell . (31b) 

Notice that /W depends only on the internal mass m 2 , whereas f( 2 \m 2 ,m 2 ) depends on the 
external mass m 2 as well. 

Starting from the definition of Aa k in (12), we have, 

(SZ ak - SQ ak ) a k = Aa k - bZ^, k a k - 8Qa k - ^8Za k , (32) 

where all the quantities on the r.h.s. have already been computed. The terms in SZ^ k and 8Q 
on the r.h.s. already have the form required to determine the l.h.s.. In order to cast the other 
ones as a dilatation times a rotation, we proceed as in appendix D. Thus, for the wave-function 
renormalization matrix we write, 

SZa k = (a k -8Za k ) a k + ((8Za k ) A a k ) a k . (33) 

with the usual definition of the wedge product as uAv = u0v — v0u. For the vertex counterterm 
we have (Aa fe ) c = <5ia fc (a k ) c + (8 2 a k ) c with, 

2 N s Ns 

5ia k = 0^ " E K)e 2 / (1 V?), (S 2 a k ) c = -(a k ) c ^2(a k ) d 2 fW(mlm 2 ). (34) 

c=l d=l 

Defining the diagonal matrix (5B) a b = —5 a b J2d=i ( a k)d 2 f^H'm 2 , n^d) we can wr it e 8 2 a k as 8 2 a k = 
(a k -8Ba k ) a k + ({8Ba k ) A a k ) a k , whence, 

Aa k = (Sia k + a k -8Ba k ) a k + ((8Ba k ) A a k ) a k . (35) 

Notice that only 5\a k is divergent. Gathering all contributions together we get, 

8Z ak = 5 1 a k + {a k -8Ba k )-5Z i , k -]^a k -8Za k , 8Q ak = 8Q - {8Ba k ) Aa k + ^(8Za k ) Aa k . (36) 

These constants determine the renormalization of a k and a k separately, the latter being necessary 
for the renormalization of V4. 

2.4 MS scheme 

In MS scheme we set 8Q = by definition. We choose our flavor basis so that at tree level the 
renormalized mass matrix M 2 is diagonal. This choice makes the tree-level propagators diagonal 
which, although by no means mandatory, conveniently simplifies the treatment. We write M 2 = 
M' 2 + M 2 , with M' 2 = dmg(m\, . . . ,m 2 Ns ) containing the renormalized masses, and M 2 the 
off-diagonal elements, which are of second order in the coupling constants, and are treated as 
interaction terms. As pointed out in [1], since the mass matrix is symmetric, we can always 
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write M 2 = e-x.p(E)M' 2 exp(— E), with M' 2 the diagonal matrix of eigenvalues and exp(2£) 
the orthogonal matrix of normalized eigenvectors of M 2 . E is antisymmetric and, due to our 
choice of tree-level flavor basis, it is of second order in the coupling constants. At one-loop level, 

M 2 = M' 2 + [E, M' 2 \ H , with M 2 = [E, M' 2 \ having vanishing diagonal entries. 

We rewrite the cf> field Lagrangian as, 



with, 



= -^>-(n + M' 2 )(f> - l -4>-M 2 4> - X -<S>-{8ZU + AM 2 )0 (37) 



AM 2 = SM 2 + [SQ m , M 2 ] + ^{SZ, M 2 }. (38) 



Counterterms to the <f> two-point function are obtained from (17) as, 

N f 
k=l 

8Z is invariant under — > — o^, fe = 1, . . . , iVj, as expected and, in the notation of (5), <5Z^fc = 
-a|/(47r 2 e), <5^ fc = = 8Z' . With these values (39) for SZ and AM 2 we can find the mass 
renormalization constants 8M 2 and SQ m from (12), 

i 2,2^ 

*M 2 6 + [SQ m ,M\ b = (/f 2 ) afc + (n£) + ^EKW«^- ( 4 °) 

fc=i 

Unlike the case in OS scheme, in MS SQ m is not needed in order to renormalize the coupling 
constants. We use below the complete mass matrix renormalization as given by (40) to derive 
the one- loop RG equation for M 2 . Regarding the choice of M 2 and its inherent ambiguities, 
considerations completely analogous to those given in Section 2.4 of [1] are valid, there is no need 
to repeat them here. 

The cj) self-energy in this scheme is then given by the e- independent terms in (17), augmented 
by the off-diagonal mass matrix, H(p 2 ) = (^(p 2 )) finite + M 2 . The relation among renormalized 
masses and pole parameters is obtained from T2 aa {s Pa ) = 0, with s Pa = m 2 a —im Pa T Pa . At one-loop 
level, 

N t 

m 2 a = m 2 pa - Xtn aa (m 2 pa ), T Pa = ^ £ (a k ) a 2 . (41) 

k=i 

The one-loop self-energies and wave-function and mass counterterms for \ an d fermion fields are 
given by, 

1 1 Ns 

6Z ^=-jQ^; a li n^ fc (p 2 ) = 32 oa Yl ( ak ^ 2 ( m c Q o( m c) + (P 2 ~ m 2 c )b (p 2 , m 2 , 0)) , (42) 

c=l 

1 1 Na 

5Z X = 0, bm\ = — 2~ Tr (if 2 ) , U x (p 2 ) = £ H ab H ba b (p 2 , m 2 , m 2 ), (43) 



a,b=l 



from whence m 2 = m 2 p — n x (m 2 p ). This expression satisfies the RG equation for m 2 obtained 
from its counterterm, ~pdm x /dj2 = l/(167r 2 )Tr (ff 2 ), as can be checked from (43). 
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The unrenormalized 4>~X vertex function has the form given in (22). The only divergent con- 
tributions to r a fr(pi,p2) at one loop come from diagrams with one V4 vertex. The detailed form 
of A.H in the case N s = 2 is given in appendix A. The renormalization constants 8Zh and SV 
defined in (7) satisfy (26) with 8Q = 0. If we set all a k = and V4 flavor-isotropic we can set 
SV = 0, since in that case SV can only depend on H and g and therefore commutes with H. In 
the case a k 7^ we consider in this paper, SV 7^ must be determined from (26). 

We consider, finally, the renormalization of Yukawa couplings. Since SQ = and SZ = 
Ylf=i SZj^&j <g) &j, we can write the counterterm to the Yukawa interaction as, 

N f 

Aa k = Uz ak + 6Z fk + -Szf j a k - 8Q ak a k + - ^ Szf (a r a k ) a r (44) 

In order for the interaction to retain its form under renormalization, we must have SQ ak a k € 
S Y - This is guaranteed in MS by the softly broken 1j2y symmetry. From the expression for the 
unrenormalized Yukawa vertex scalar form factor given above, we obtain, A.a k = 1/ (87r 2 e) a\a k . 
Substituting this value and the values from SZ^ k and 5Z^ in (44), 

N f 

5 2 a k 



(5Z ak - SQ ak ) a k = Y^4«fc + ^T e zZ ^jk&j, 



3 1 Nj 1 ( Nj \ 

3=1 \i =1 / 



(45) 



where Xj k = &j-a k and the second line is obtained by substituting aj in the first one by aj = 
x k j& k + (aj A a k ) a k . Clearly, 8Q ak is determined only up to an antisymmetric matrix nullifying 
a k . From (45) we see that we can set, 

N f N f 

3 1 \ 1 \ 1 

SZak = 16^°' + 8^5Z a ^ fe ' SQ ak = --^Y, a2 J x ^( & 3^^k) = ^(SZa k )Aa k , (46) 

i=i j=i 

which fixes the renormalization of the Yukawa couplings {a k } and versors {a k }. If the {a k } are 
mutually orthogonal, SQ ak = 0, and SZ ak = 5a|/(167r 2 e) does not depend on aj with j / k. In this 
case the fermion sector of the model reduces to Nf copies of the case Nf = 1. Similar simplifications 
occur if all {a k } are collinear, or if they are divided into mutually orthogonal subsets of collinear 
versors. We discuss these geometric issues in section 2.5, in connection with RG equations. The 
renormalization of the potential terms V^{(j)) is discussed in appendix A. Notice that only the 
matrices SQ ak enter the counterterms to V4, but not the constants SZ ak . 

The RG equations for the Yukawa couplings {a k } and versors {a k } have the form, 



-da k _da k 

^=Pa k , M^=/V (47) 



From (8), with the renormalization constants (46) and retaining terms through 0(afy, we get, 

N f 
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We see that if the coupling vectors are pairwise orthogonal, Xjk = if j ^ k, the RG equations 
for the ak decouple. In the general case the evolution is coupled, but it does not depend on the 
sign of Xjk- In particular, (3 ak > for all k at d = 4 so all Yukawa couplings grow with log(Tl) as 
expected. For the Yukawa versors we obtain, 

3=1 

We discuss this equation in more detail below. 

The RG evolution of the coupling constants in V4 is needed in order to obtain the RG equations 
for H and M 2 , but only to lowest order Jldg/dp = —eg + 0(g 3 ), and analogously for g^j and gijki- 
From (7) we obtain a RG equation for H with the renormalization constants SZh, SV as given 
by (26) with SQ = 0, 

-p^r=(3 H , P H = -^H + e f i-^ 2 AH+-± r J2{a k ®a k ,H}. (50) 
^ 71 e k=i 

AH gathers the divergent part of the contributions from V4 to the <j>~x vertex. For N s = 2 it is 
given in (A.l). 

From (40) we find, 

Nf 

^ fc=l 

where 7 m , which has mass dimension 1, is defined by this equation. The RG equations (50) and 
(51) are algebraically complicated, due in part to the contributions to (3h and 7 m from diagrams 
with one insertion of V4. We will only make the following simple remarks about these equations. 
If {a.^} are eigenvectors of M 2 at a given scale, they must lie on its orthogonal eigenspaces which, 
since M 2 is regular, must be one-dimensional. Therefore, any two versors aj, and dj must be either 
collinear or orthogonal. As discussed below, such configurations are fixed points of the RG flow. 
Thus, if we require M 2 au = m? a* for all k at all scales, differentiating we obtain j m ak = ^a k &k 
for some eigenvalue Xa k , leading to [iW 2 ^^] = in S Y - In short, {a^} can be eigenvectors of M 2 
at all scales only if M 2 and 7 m can be diagonalized simultaneously in S Y for all ~p. The converse is 
clearly valid. If some m 2 are degenerate, the versors can be eigenvectors of M 2 without being 
at a RG fixed point. Analogous considerations hold for H and (3h- 

The anomalous dimension of x vanishes. For the fermion fields, with 8Z^ k from (42) we find 
7V>* = — V(327r 2 )afc- From (39) we get for the anomalous dimension of <j>, 

N f 

1 d 1 

k=l 

7 is symmetric, which is consistent with SQ = 0. Notice that at this order (f> components in 
directions orthogonal to S Y have vanishing anomalous dimension, and that only j3 ak contributes 
to 7, but not /3 &k . This prompts us to look at the anomalous dimension of • 0, which receives 
contributions from both /3 a and 7. Those two contributions partially cancel each other and we 
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get, symbolically, 

d 1 Nf 

~p-j=a k -4> = r] k -4>, rj k = ^a)x 2 jk a k . (53) 

j'=i 

At this order the RG evolution of the components a k ■ 4> of <j> is decoupled from the evolution of 
with j / A;. 

2.5 RG flow in flavor space 

The RG equations for the Yukawa coupling constants {a k } as given by (47) and (48) show that 
those couplings grow monothonically with log (/I). Since (3 ak > for all k, there can be no fixed 
points. In this section we consider the qualitative analysis of the RG flow of the Yukawa coupling 
versors {a k }. Specifically, we discuss its fixed points and invariant manifolds, and their stability. 

We find it convenient to describe the problem not in terms of the {a k }, but of its associated 
Gram matrix x G R N f xN f, x,- L j = ai &j. x is symmetric, with diagonal entries equal to 1, and 
rank(a;) = d Y - We consider Xij with 1 < i < j < Nf as our basic variables. If the quantity Xji with 
i < j occurs in an equation, it is to be interpreted as xji = x,^ both in its value and as a variable. 
In the same sense, /3 Xji = (3 Xij if i < j, and xu = 1, (3 Xii = 0. There are then ( 2 f ) independent 
variables. 

Prom (49) we have, 

Nt 

7^=^ = Ac, f3 Xjk = -^Y^af {%ik{xij - x ik x kj ) + xij(xi k - xijXjk)) • (54) 

By definition, \xij | < 1 for all If Xij — ±1 for some i,j, from (54) we see that f3 Xij vanishes. 
Thus, (54) is consistent with \xij\ < 1. 

We are interested in the fixed points of (54), zeros of f3 x , which are independent of the values 
of ai > 0, I = 1, . . . , Nf. Those are given by all (^-tuples x° = (x°ij) with i < j, such that all 
x°ij = 0, ±1, and satisfy the geometric consistency condition, 

if x°ij = ±1 then x^jx"^ = x° ik for all 1 < k < N f . (55) 

The meaning of this condition is apparent: given that &i and a,j are parallel, if a k is parallel (resp. 
antiparallel, orthogonal) to aj, then it must be parallel (resp. antiparallel, orthogonal) to dj, and 
analogously if dj and a,j are antiparallel. 

If x° is a fixed point with some components equal to -1, then x ' with x '^ = \x°ij\ is also a 
fixed point, as can be seen from (55). Since we can always flip any a k in C by reparametrizing the 
corresponding fermion field as if) k — > J^ipki x ° an d x °' are equivalent. In particular, the linearized 
(3 X function has the same eigenvalues and eigenvectors at x° and x ', so we need consider only 
those fixed points x° with x°ij = or 1 and satisfying (55). 

The fixed points can be grouped in three types. First, there is the fixed point x° with x°ij = 
for all i < j, corresponding to the Yukawa versors {d&} being pairwise orthogonal. Clearly, in this 
case we must have d Y = rank(a: ) = Nf < N s . The linearized (3 X function is diagonal and strictly 
positive definite, 

8n 2 d _pl = y ( 2 + 2) > . (56) 
This fixed point is therefore IR stable. 
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The second type comprises fixed points with Xij = ±1 for all i < j, i.e., all collinear, d Y = 1. 
These are all equivalent, in the sense defined above, to x° with x°ij = 1 for all i < j. The linearized 
(3 X function is scalar and strictly negative definite at these fixed points, 

N f 

= -2b ik b j{ V a 2 m < 0, (57) 
— i 

m=l 

which are then UV stable. 

The third type of fixed points consists of those x° with not all x ^ = and not all x°ij = ±1. In 
this case 1 < d Y = rank(a;°) < Nf, and only those x° with rank(a; ) < N s are possible. These fixed 
points correspond to the {&k} lying along r orthogonal directions in flavor space, with 1 < r < Nf, 
so that at least two of them are collinear, but some of them are orthogonal to each other. As shown 
in appendix E, these are all saddle points. 

The set of fixed points just described are precisely the points in parameter space where C is 
approximately invariant under the softly-broken discrete symmetries Z^y, etc., discussed at the 
beginning of this section. An important particular case occurs when at a fixed point of second 
or third type not only some coupling vectors are collinear, but also their magnitudes are equal. 
For example, if otj = aj, i ^ j, and they are orthogonal to all other coupling vectors, then 
from (47)— (49) we see that the equality is maintained at all scales. In that case, C is invariant 
under an exact SU(2) symmetry acting on the fermion fields ipi, ipj, on top of the corresponding, 
generally approximate, discrete symmetry. If there are 1 < r < Nf groups of i%i > 1 equal coupling 
vectors, with each group orthogonal to all other coupling vectors, the fermion sector has an exact 
SU(rii) x • • • x SU(n r ) invariance. 

Going back to the general case, not all of the fixed points described above can be reached 
for given values of Nf and N s , and for given initial conditions for (47) or (54). For instance, as 
mentioned above, the fixed point of first type is not reachable if Nf > N s . More generally, we are 
led to ask how linear dependence and independence of {&k}, and d Y , evolve under the RG flow. 

The Yukawa versors are linearly independent if and only if x is non-singular. From (54), taking 
into account Xij = Xji and xu = 1, we obtain the RG equation for det(a;), 

Nt 

d 1 

JL— det(aj) = 7 d det(a;), j d = --^ ^ a\x\ < 0. (58) 

We see that if at a scale 7i we have det(a?) = (resp. > 0, < 0), then det(a:) = (resp. > 0, 
< 0) at all finite scales 71. Thus, linear (in)dependence of {a^} is preserved by the RG flow (but 
not asymptotically preserved, however, since at the UV fixed points d Y = 1). Assume that at 
7i the set {a^} is linearly independent. We denote by C Y the "Yukawa simplex" formed by all 
convex combinations of {&k}, whose volume is proportional to det(aj). From (58), det(a;)|^ = 
(7i/7i ) 7d det(a;)|jj at one loop, so that C Y loses volume as 71 grows at a rate dictated by the 
anomalous dimension 7^ of det(ai). We notice also that the r.h.s. of (58) vanishes at all fixed points 
because det(cc) vanishes in all cases, except for the first-type fixed point, at which 7^ = 0. 

We can, further, consider the cofactors of order k, j N k (x; Nf,k), < k < Nf, 1 < 

ii, . . . ,iN f -k,ji, ■ ■ ■ ,jN f -k < Aj, defined up to a sign as the determinant of the k x k matrix 

obtained from x by deleting rows %\ iN-k and columns j\, ... ,]N-k- Cj i Jjv _ k (x; Nf,k) can 
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be written as 

^ N N 

^h-JN-k ^) = j^l £ m 1 ...m k i 1 ...i N _ k £n 1 ...n k j 1 ...j N _ k X mi n 1 ■ ■ ■ X mk n k - (59) 

' mi,...,m k =l ni,...,n k =l 

Since a; is symmetric, in particular diagonalizable, we are interested only in diagonal cofactors such 

_ k 

as j _, (a;; Nf, k), which is the determinant of the Gram matrix associated to {&k}k^h,...,i N _ fc - 
In particular, C(x;Nf,Nf) = det(a;). For a diagonalizable matrix such as x it is easily seen that 
rank(a?) = d Y (with <i y < Nf) if and only if iN _ k {x; Nf, k) = for all k > d Y and for all ii, 

. . . ,Z7v f _fc, but C- ■ (x:Nf.d Y ) ^ for some set of indices i\ . . .%N f -dv 

Using (54) and (59) and taking into account that and xu = 1, we obtain, 

N f -k ^ 

Ti-lnC. ■ f \N f ,k)=-f, f a . f "(N f ,k)--^ Y a 2 t C Z Nf (Nt,k + 1), 

^ m=l 1 



m=l 



V x 2 - 1 



n=l 

\n^U,— ,iN f -k j 



where the argument x has been suppressed for brevity, and the caret on an index indicates that 
index is to be omitted. Setting k = Nf in (60) we recover (58). 

Assume now that rank(a:) = d Y < Nf at some scale /Z . Then C(x; Nf, Nf) = at all scales by 
(58) and, from (60), the (Nf — l)-cofactor satisfies, 

ji±Cf(x;N f ,N f - 1) = 7 % f ^Cl(x;N f ,Nf - 1). (61) 

If d Y = Nf — 1, then at /x we must have Cj(x; Nf, Nf — 1) ^ for some j. From (61) we see that 

in that case Cj (x; N f , N f - 1) ^ at all finite scales ~p. If d Y < N f - 1, then Cj (x; N f , N f - 1) = 
at ~p for all j, and by (61) they vanish at all finite scales. Iterating this argument, we conclude 

that if at some scale ~fi we have C {i i N _ k {x; Nf, k) = for all ii, . . . ,iN f -k an d all k > d Y , and 

ijy d (a;; iVy, d y ) 7^ for some set of indices i±, . . . ,ij^ f -d Y , then that situation persists at all 

finite scales Ji. Therefore, rank(a?) = d Y is RG invariant (though not asymptotically invariant). 

For k < Nf, 7 fc 3 can be positive, negative or zero, depending on the values of and 
Xij. In a small neighborhood of an UV fixed point, x^ will be close enough to 1 for all so 

il...i^ f _ k 

that 7^ 3 < for 1 < k < Nf. In that situation, (60) with k = d Y describes how C Y loses 
ciy-dimensional volume as Ji grows. Needless to say, we are implicitly assuming that even though 
x is near an UV fixed point, possibly as a result of a choice of initial conditions, the couplings ak 
are small enough for perturbation theory to be valid. 

3 An S*7(iV)-symmetric model with SSB 

In this section we turn to an SU(N) symmetric model with scalar fields in the adjoint represen- 
tation, <p = (<pi, . . . , <j>jj), and fermion fields in the fundamental representation, i/> = (ipi, . . . , tpx)- 
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We denote n = n 2 — 1 for brevity, A = (Ai, . . . , A^y) for the vector of Gell-Mann matrices as defined 

— * — * 

in appendix F, and cf) = The classical Lagrangian in the broken symmetry phase of the model 
has the form, 

C = -lf{D- u 2 )$- V A {$,N)+^i M - g$<W, (62) 

with v 2 > and 14(0, N) the most general SU(N) invariant quartic homogeneous polynomial in 
cf), as given below. C is invariant under the U(l) symmetry related to fermion number conservation, 
the %2 chiral symmetry cf) — > —(f), if) — > 75^, and the SU(N) transformations, 

N _ 

$ _> e -^-A/2 V , ; a ^ ^ iU(#K i?a 6 (0) = (e^ X/2 A a e-^ X / 2 A b ) . (63) 

6=1 

The SU (N) currents are given by j£ = tp^K/^ip-J2ab fcabd^fiaA, with f cab the SU (TV) structure 
constants. Both Z2 and SU(N) symmetries are broken if (cf)} / 0. 

We can build two quartic SU(N) singlets out of cf), (Y^ij (pijPji) 2 and Y^ijki <f>ij<f>jk<f>kl<i>li, giving 
rise to two terms in V4. Given two vectors a, (3 in the adjoint representation, we introduce the 
notation (a V/3) a = Ylbc=i d a b c ab(3 c , with <i a b c the anticommutator constants for A matrices. With 
these notations, V4 can be variously written as, 



Vi(<p.X) - '- (f) - '— ( (p / cf> ) 7 I -j dabhd cdh (f) a (f) b (f) c cp d 

abcdh=l (64) 




TV). 



The coupling constants A, A' are defined to be positive. The choice of sign for the term in A' in 
(64) will be explained shortly. It is easily seen, however, that for appropriate ranges of positive 
A and A', V4 is bounded below. Given any Hermitian matrix A, not necessarily traceless, we can 
evaluate V±{A, N) in a flavor basis such that A = diag(«i, . . . , ajv), with aj real, 



via *n 1 X JV-IA^ 4 1 /A A'\ A 22 

VAA,N) = -[ A > a* + - -H > a-a 2 . 

v ' ; 8V2 /V J 8V2 N ^ 1 3 

' j=l v / i^=l 



We see that if we choose A'/A < N/(2(N - 1)), then F 4 (A, TV) > for all Hermitian A. Setting 
A'/A < 1/2 we ensure that V4 is bounded below, and in fact positive, for all N > 2. 

-* -* -»2 

We denote the local minima of the tree-level potential V(cf), N) = V A (cf>, N) — v 2 /2 cf) by Vt \ = 

#(o)-A = U(o)V(o)-A. Due to ££/ (AT) invariance the manifold of minima of V(cf>, N) consists of a disjoint 
union of SU(N) orbits, each of which contains several diagonal matrices with the same eigenvalues, 
related to each other by SU(N) transformations which are reorderings of the flavor basis. Each 
one of these orbits constitutes an SU(N) equivalence class of minima of V, or "modulus" for short. 
We can choose a diagonal matrix v',^ = w' (0) X^L 2 (*(o) witn ^""(o) = J2n=2 (V (0) )l = 1, out 
of each orbit as representative of the corresponding equivalence class. Thus, we can parametrize 
classical moduli space as a submanifold of the algebra of diagonal matrices of su(N) generated by 
{\n},n = 2,...,N. 

For N = 2, d abc = and V(cf>, 2) is isotropic in flavor space. Its minima are all vectors V( ) with 
v (o) = vj\f\, all of them lying in the SU(2) orbit of v'^ = v^X^. The space of moduli consists 
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in this case of a single point. There are two Goldstone bosons and a massive scalar with tree-level 
mass \^2u, and the two fermion fields are mass-degenerate. 

The case N = 3 is special. Using the fact that is a root of its characteristic polynomial and 
that Tr(» = 0, we find that Tr(0 4 ) = l/2(Tr(0 2 )) 2 [4]. In this case also V A {$,3) is isotropic, 
V 4 ((j>,3) = 1/16(A - A'/3)(Tr {4> 2 )) 2 ■ (We notice, incidentally, that we can take A'/A as large as 3 
without losing positivity of V^{4>, 3).) From this expression, and (62), we see that at tree level the 
scalar sector of the model is invariant under the action of the fundamental representation of the 
group SO(8). The classical minima of the potential are the vectors t>( ) with = vj \J\ — A'/3. 
Once a direction for vr \ is chosen its stationary group is 50(7), so at tree level we have seven 
massless and one massive scalar field with mass \[2v. Clearly, 50(8) is not a symmetry of the 
model, since it is explicitly broken by the fermion sector. In particular vevs V(q\ with different 
eigenvalues, which correspond to different moduli, lead to different fermion mass-spectra at tree 
level. 

Moduli space can be parametrized by the circle {v', J9) = v^(cos9X^+sm9Xs)}. Since v', J9) 
and t?/QN (9 + 27r/3) differ only by a permutation of their diagonal elements, if we want to factor the 
action of the permutation subgroup Z3 of SU (3) we must restrict the range of 9 to — ir/3 < 9 < ir/3. 
All V(o)(0) with 9 = tt/2 + mr/3, in particular t>'^(7r/2) = As, have two equal eigenvalues. Thus, 
their stationary subgroups are isomorphic to SU(2) x U(l). In this case there are four Goldstone 
bosons, three scalar fields acquire mass radiatively, and tree-level fermion mass ratios are 1:1:2. 
On the other hand, v'^(9) with 9 7^ tt/2 + n-7r/3 have three different eigenvalues, and a stationary 
subgroup isomorphic to U(l) x U(l). There are six Goldstone bosons in this case, with one scalar 
acquiring mass radiatively. All vacua are degenerate at tree level. 

For N > 3, V^(cf),N) is not isotropic in flavor space. Like in the previous case, we can 
parametrize moduli space by a set of minima v', -. of V4 within the subalgebra of diagonal ma- 
trices of su(N). Unlike the case N = 3, for N > 3 moduli space consists of a finite set of isolated 
points. Rather than trying to map all of moduli space for all N > 3, in this paper we restrict 
ourselves to pointing out that for all N > 2 there is a local minimum of the potential of the form 
Vr \ = v^Xjj. This minimum has the largest possible stationary group, SU(N — 1) x U(l), leading 
to the fewest possible Goldstone bosons, 2N — 2. 

Evaluating the extremum conditions dV(<f>, N)/d(j)t = 0, t = 1, . . . , N, for the tree-level poten- 
tial at f (0) = V( 0) \ w (i.e., («( ))a = u (o) we obt ain 



"^Tx^w Dsd ^ = - iN - 2) \jNTN^Ty (65) 

where use was made of the properties of the constants d a b c with one or more indices equal to N, as 
given in appendix F. With the assumption A' < A/2, the argument of the square root is positive 
for N > 3. At this extremum second derivatives d 2 V /d4> s d<pt with s ^ t vanish. For s = twe find, 



d 2 v 



o d 2 V 
= 2v z > 0, and — T 

<t>= v (0) x W 



N — 3 

= 2X ' v %)^-r T if 2<n<iV. (66) 



(0) N _ 1 

<£=f(0) *]v 



The r.h.s. in the second eq. in (66) vanishes if N = 3, and it is positive for N > 3 as long as 
A' > 0, which explains our choice of sign in (64). Due to SU(N) invariance, there is no need to 
evaluate the second derivatives in the directions of the non-diagonal generators: d 2 V/d(p 2 with 
1 < t < N — 1 — 1 and t 7^ n leads to the same result as in (66) with t = n, n < N. If we 
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set N — 1 < t < N, we get d 2 V/d(f> 2 = at vt y as required by Goldstone's theorem. Thus, 
v (o) = v (o)Xjt 1S a minimum of the tree-level potential for N > 3 if A' > 0. If A' < and N > 3 it 
is a saddle point. We may also ask whether there are other minima oc Xjt, with 2 < n < N. This 
occurs only for N = 3. For N > 3, there is an extremum of the potential only in the direction of 
A3 but numerically it turns out to be a saddle point for N = 4 and 5. All other minima in the 
subalgebra of diagonal matrices of su(N) must be linear combinations of more than one A^. 

Shifting the field = ip + v^y and substituting in the classical Lagrangian (62), we obtain the 
tree-level masses for the scalar fields 



at _ 3 

m\ = 2X'vf 0) — — - if 1 < a < N - 1, m\ = if (N - l) 2 < a < N, m^= 2v 2 . (67) 



For N = 3, as expected, m 2 = for a < 7, (£1,2,3 acquiring mass at one loop. For the fermion 
sector we have the mass matrix = gv^Xjj or, more explicitly, 



m^=m^ if 1 < i < iV — 1, m^, N = —(N — l)m^, = 9V(o)^ j^jj _ jj - (68) 

Notice that, at tree level, X)£Li m V>, = an< ^ S^Li m % i = 2g 2 v 2 Q y For the purpose of the pertur- 
bative expansion we consider A' to be of the same order of magnitude as A, and g = ©(A 1 / 2 ), so 
that = O(l). 

In the quantum Lagrangian we introduce renormalization constants for the parameters, v 2 = 
v 2 + 5u 2 , A = /J, e (X + 5X), X' = //(A' + 5X') and g = ii e/2 Z g g. We restrict ourselves to MS 
scheme. Thus, wave function renormalization matrices reduce to the scalars 5Z for cf> and SZ^ for 

due to SU(N) symmetry. We write C in terms of shifted fields (f> = <p + /U~ e / 2 0, where fi is the 
MS mass scale, e = 4 — d, and the vev is given by = 0( O ) + vn\ with 0( O ) the tree-level vev, 
V(oyX = v^Xjj, and 0^) = 0(X 1 ^ 2 ) to be determined by the requirement that (0|<p|0) = at one 
loop. 

With these definitions, and retaining only terms contributing at one-loop level, C is written as, 

C = C f + C 1 + C 2 + C hlt + C ct (69) 
Cf = —-(p-Dtp + \p-M<p + ifti flip — tpM^ip 

d = /i _ 3 U 2 V {l y0- /i"3 Xvf 0) V (1 y<p - 2/i~5 A (#(()) ■#(!)) (0( O )-<£) 

+ A' (0( O) V (o) ) • (0(1) V <£) + 2// _ i A' (0 (o) V 0(i)) • (0(0) V <p) 
C 2 = -X (0(0) -0(i)) 2 - 2A (0(o)-<?) (#(i)-£) + X' (0(o) V (1) )- (<pV<p) 
+ 2A' (0(o) V <p) • (0 (1) V <£) - gipv {1 yXip 



Ant = — A*2 A (0 (o) - 0) (f + /i2 A' (0(o) V <£) ■(<£ V (p) - fi 6 ^ (ip 2 ) 2 + fi 6 ^ (<pV Cp) 2 - g^tpip-Xil) 



1 _ - e 

£ c t = — -8Z(p-U(p + SZ^xfri flip — gAZ g ipV( yXtp — gfi^ AZgipip-Xtp 



+ H 5 ( A^ 2 - t; 2 0) A A) 0(o) • £ + /i 3 A A' (0(o) V (o) ) • (0 (o) V 0) + - ( A^ 2 - t; 2 0) A a) <^ 
AA' 

- A ^ (0(o) -(p) 2 + — (0(o) V (O ))-(^V cp) + AA' (0(o) V ip) 2 - fi^AX (v (0 y<p) Cp 2 
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+ fii AX' (v {0) (<ff + fi e ^- {$V0? , 



where is given by (68), (v( )) a = v^S al7 with v {0) from (65), (p-M(p = Yf a =\ m lfl with m l 
from (67), and 

AA = 5X + 2X5Z, AX' = 6X' + 2X'8Z, Av 2 = Sv 2 + v 2 5Z, AZ g = 5Z g + 5Z^ + hz. (70) 

The Lagrangian C x = C(A 1/2 ), and C 2 = O(X). 

The condition that the finite part of (0|<^(0)|0) must vanish at one loop leads to (v^) a = for 
a < N and, 

^ 2 (Z(i))w = -I^™?°o K) ^(o)iV(iV - 2) (a + X %* N Z*l ) ~ J^ V{0)m ^ a ° ( A 
x/ 2(A-2) 2 \ A( 0) 1 , 2 . A( 0) (A-1) 2 , 2 , 

iv(iv^i) J + -^iv^T) °o K) + ^Sv^ ao KJ • 

At A = 3, depends on A and A' only through A — A'/3, as it should. The requirement that the 
divergent piece of the ip 1-point function vanishes yields the relation 

a 2 2 ax 2 2(iV-2) 2 , 3 2 / 2(A-2) 2 \ 

A, 2 - V 2 0) AX + ^oj^H) AV = (A - ^^yA'j 

1 , (A-3)(A-2)A ( 2(N - 4) A 2 A 2 - 3 A + 3 2 4 ^ 
" A^l A [ X + N(N-1) X ) + Wt A(A-l) » 5 • (?2) 

Together with two other similar equations from the 2-point function, (72) determines the countert- 
erms on the l.h.s.. 

The self-energy function for the (p field must be symmetric in flavor space, due to CPT invari- 
ance, and invariant under the SU(N — 1) x U(l) symmetry of C and v. Those two requirements 
together imply that it must be diagonal, of the form diag(a, . . . , a, b, . . . , b, c) where the first block 
is A — lxJV- 1, the second (2 A — 2) x (2 A — 2) and the third lxl. More generally, the triviality 
of mixing in this model holds for any pattern of symmetry breaking. Any VEV of the scalar field 
has a stationary group containing a subgroup of SU(N) isomorphic to {/(l)^ -1 . This symmetry, 
together with the fact that the scalar field self-energy must be symmetric in its flavor indices due 
to CPT invariance, is enough to make the self-energy function diagonal in the basis in which the 
VEV is. 

At one loop the self-energy function can be written as, 



iW 2 ) = n<JV) + nilV) + £n£-V), (73) 

nSV) = Sab (2Xv {oy v {1) - Av 2 + AXv 2 0) ) + 8^2 (2Xv (oy v {1) + AAtf 0) ) 
- 2 (v (Q yv {1) X' + ^AX'] (dj^ + 2<% J d^Jab - P 2 S ab SZ, 



,2 

2 

N \ N 



U ab (P 2 ) = 6 abX [ Yl A ^ m2 d) + 2 Mm 2 a) \ ~ ( d lnndlab + 2d lna d lnb ) A)«), 

,d=l J l,n=l 
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nifV) = 2^ )A 2 { 5 ab 2B (p 2 ,m 2 w ,m 2 a )+5 a - N 5 b - N ( 2B (p 2 ,m2,m2)+ 



, c=l 



+4B (p 2 ,mL, mL^ 



at 



ifi'V) = 2^ 2 0) A' 2 ^ {% aa % &6 + 2d Waa d md + 2% cc d m6 + 2(d : 
c,d=l 

+ M jv cc <%<ta} d a cddbcdBo(p 2 ,ml, m 2 d ), 



' N N 



^ab^iP 2 ) = - 2v fo) XX ' { 6 aN^2 d bcc{djj bb + 2d Jrcc )B (p 2 , m 2 ., ml) + 5 b lv^da CC ((%aa 

c=l c=l 

V 

+2d 7rcc )B (p 2 ,m 2 c ,m 2 c ) + 4c% a6 (c% aa + <%» + <%a7]v) A)(;P 2 , rn^, m 2 a ) j , 

AT 

n £" 4 V) = -V A^Aft.. (A)("4) + + ((m^ + m v ,.) 2 -p 2 ) B (p 2 ,mJ.,mJ.)) . 

The sums over flavor indices can be computed explicitly with the aid of the formulas in appendix F, 
the resulting self-energy Tl ab (p 2 ) having the diagonal form described above. We skip the details here, 
to focus only on the divergent part (Ji a b{p 2 ))(x> of H a b(p 2 ) and some of the finite mass corrections. 
Requiring (H aa (p 2 ))oo = for the three ranges 1 < a < N — 1, N — 1 < a < N and a = N, leads 
to three linear relations among Ai/ 2 , A A and AA', one of which is identical to (72). After solving 
those relations we find, 



(74) 

^ _ O + 1)A _ 2 ^ZM , A y _ 1 (-4^A° + 12AA<) . 



For N = 3 these values for AA and AA' are arbitrary, since the relation (n aa (p 2 ))oo = depends 
only on AA — AA'/3. It is only this last quantity that is univocally determined and it, as well as 
Az/ 2 , depends only on A — A'/3. 

The finite mass corrections can be obtained explicitly from (73). There are two particular cases 
of interest. One is the calculation of II aa (p 2 = 0) with N — 1 < a < N whose vanishing, required 
by Goldstone's theorem, constitutes a non-trivial verification of the result (73) for the self-energy. 
We omit the details of that calculation for brevity. The other interesting case is the computation 
of U aa (p 2 = 0) with a = 1, 2, 3, for N = 3, since that is the mass for 921,2,3, which are massless at 
tree level. We again omit the details and only quote the result, 

n Qa (p 2 = 0) = J^f - (l log(2) - I) , a = 1,2, 3, N = 3. (75) 



7r 2 (3A — A') V3 

This is the common value of the renormalized masses w- 2 2 3 , which at this order are also pole 
masses. If g = those masses vanish, as expected, since the Lagrangian without fermions has an 
50(8) symmetry leading to seven Goldstone bosons. 
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For the fermion self-energy we have, 

Uij(p) = U%\p 2 ) p 1 + Iljf (p 2 ), l<i,j<N, (76a) 

N N 



f 1 1 1 

n — 1 1 — 1 V 



32vrV 



a=l Z=l 

x (m^m 2 ) - m^a (m^) + (m^ - m 2 )6 (p 2 , m 2 ^ , m^)) }> , (76b) 

Af AT 



njf (p 2 ) = gAZ g (v {oy \) . . + <? -X) . . + ^ X;(Aa)«(Ao)« 

U %3 a=l J=l 



16vr 2 



fco^m^m 2 )} , (76c) 



from whence, 



q 2 N 2 -l q 2 

SZ * = -A— AZ » = iAv? (77) 

From the renormalized self-energy we can find the relation between the renormalized and pole 
fermion masses. 

In order to complete the renormalization in MS we have to consider the RG evolution of the 
Lagrangian parameters and fields. From the renormalization constants given above, we find, 



W--eA-— g A + ~qA <7 + r— 2 AA , (i g - --g + 



(78) 



2vr 2 TV ' vr 2 y ' 2vr 2 ' ^ 2* ' 16vr 2 N 3 ' 

At d = 4, /3 9 > and, for TV > 3 and A'/A < 1/2, 

At2 + 7 1 JV 2 -4 . 

so there are no non-trivial fixed points for these coupling constants at this order. There is an 
invariant manifold, g = 0, A' = 0, A > 0, corresponding to the SO(N) invariant A(</>) 4 theory 
without fermions. For N = 3 only A — A'/3 is meaningful and from (78) we get, 



>-AA)^Aj) 2 ^AA 



which is also positive at d = 4 if A'/A < 3. The anomalous dimensions for the mass parameter v 2 
and the fields are then as follows, 

1 fN^ + l A^-4 2 \ g 2 g 2 N 2 - 1 

^ = 2^{-^r- X -^T~ X+9 J' 7 ^ = -4^' ^ = -16^^V— (?9) 

The contribution to 7^2 > from boson loops is positive for N = 2, and also for all N > 2 if 
A'/A < 1/2. 
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3.1 Explicit symmetry breaking and non trivial mixing 

A simple extension of the model (62) is obtained by substituting —v 2 in C by a symmetric N x N 
mass matrix M 2 . If we denote by P the SO(N) matrix diagonalizing M 2 , the new Lagrangian 
in the mass basis is 

C = ~$.{U + M' 2 )$-Vi{P$, N)+^i ?ij> - 5 ^A-P0V, (80) 

where M' 2 = diag(^, . . . , vj^) is the diagonalized mass matrix. P appears now as a coupling 
matrix in the Yukawa term, and also explicitly in the non-flavor-isotropic piece of the potential 
for N > 3. If P is in the SU(N) subgroup of SO(N), i.e., it is of the same form as R in (63), it 
can be removed from the Yukawa term by means of an SU(N) transformation of the fermion fields 
and, furthermore, V±(P<t>,N) = 14(0, iV), so P drops from C. That is always the case if N = 2. 
In general, for N > 3, M 2 needs not be diagonalizable by an orthogonal matrix in the SU(N) 
subgroup, and P is unremovable. In the case N = 3, V4 is SO(N) invariant and, therefore, P 
appears only in the Yukawa term. If N > 3, an SU(N) transformation of fermion fields leads to a 
coupling matrix in the Yukawa term that is different, though SU(N) equivalent, to the matrix P 
in V4. 

The Lagrangian (80) is approximately invariant under the SU(N) transformations (63) with 
R(0) substituted by P t R(0)P, only the </> mass term being non- invariant. Due to this softly 
broken SU(N) symmetry the matrix P is not renormalized in MS scheme. In other schemes 
P may receive finite renormalization, which can be parametrized as P = Q t PQ, with Q an 
orthogonal renormalization matrix (see [1] for more details on this type of parametrization) . The 
treatment in this case is completely analogous to those of the previous sections. 

4 Final remarks 

In section 2 we consider in detail a simple model with N s scalar and Nj fermion fields. The 
renormalization of the model in OS scheme is carried out in section 2.3, following the approach 
of [1]. We find the counterterms and renormalized propagators explicitly, with the exception of 
contributions from diagrams with one quartic vertex, which are given only for N s = 2 in appendix 
B. Counterterms to Yukawa couplings are also given in explicit form in (36), separately for their 
magnitudes and unit vectors. Whereas this can be done in closed form for row coupling matrices, 
it is not possible in the case of square coupling matrices. The counterterms SZh and SV to the 
Hermitian coupling matrix H are only given implicitly in equation (26). However, (26) is enough 
to renormalize the vertex and, in MS, to derive RG equations for H. Scalar fields in this model 
are unstable, so their propagators may be computed in pole-mass scheme instead of OS. We do 
not expect significant changes in the results at one loop (see, e.g., [5]). 

The renormalization in MS is considered in detail in section 2.4. Since the symmetry break- 
ing effects of lower-dimension operators do not affect wave-function and dimensionless couplings 
counterterms in this scheme, those take geometric forms as in (39) and (46). All counterterms 
and Green functions can be obtained explicitly, with the exception of H counterterms, like in the 
previous scheme. The RG equations for the Yukawa coupling vectors are analized in detail in 2.5. 
The fixed point structure of the flow is as expected from the symmetries of the model. We give 
also a stability analysis of those fixed points, and of the invariant manifolds. One conclusion of 
that analysis is the RG invariance of the dimension of the "Yukawa subspace" of flavor space. The 
results of 2.5 are independent of the masses of fermions and of other super-renormalizable couplings 
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in the Lagrangian, some of which we set to zero for simplicity. Furthermore, at one-loop level they 
are also independent of the detailed form of the quartic interactions V4, and of the term (fi-K&x 2 
discussed in section 2.2. 

There is a close analogy between the renormalization of Yukawa versors on one hand, and that of 
unitary and orthogonal matrices on the other. In both cases the couplings take values on a compact 
manifold which is covered once by the transitive action of a symmetry group, and this property 
determines the form of counterterms, as in (8) where a^o = Q ak ^k with Q ak € SO(N s ). For this 
reason we expect that the renormalization of Yukawa versors and its associated RG flow should 
be a toy model of the more complicated situations in theories with several unitary or orthogonal 
coupling matrices. 

In section 3 we consider a model with spontaneously broken SU (N) symmetry. We choose a 
vev which is a local minimum of the potential for all N, with SU(N — 1) x U(l) as stationary 
subgroup. We consider the one-loop two-point function, obtaining the masses of pseudo-Goldstone 
bosons in the case N = 3, and the MS one-loop counterterms and RG equations for all N. We 
notice that in MS the analysis of two-point functions is enough to determine the counterterms to 
the coupling constants in the quartic terms in potential, whereas in OS the opposite holds, namely, 
renormalizing the cubic and quartic terms in the potential is necessary in order to obtain the mass 
and field-strength counterterms [6]. 

Mixing of fields is trivial in this model if the symmetry is not explicitly broken, for instance 
by a mass matrix as discussed in section 3.1. The local minima of the potential strongly depend 
on the form of the mass matrix in that case. Such explicit breaking of the symmetry leads to 
the appearance in C of an orthogonal coupling matrix both in V4 and the Yukawa term. It is 
not possible to give closed expressions for two-point functions with a generic mass matrix for all 
N, as done in the case without explicit symmetry breaking. In the MS scheme counterterms to 
dimensionless couplings are independent of the mass matrix so, in particular, the coupling matrix 
P does not receive renormalization. In a scheme other than MS finite counterterms to P are 
obtained along the lines of the treatment of the previous sections and [1]. 

A The potential V 4 ((f>) 

In section 2 the quartic self-interactions are denoted V^cf)). In this appendix we give explicit 
expressions for one-loop contributions to 2- and 3- point functions from diagrams containing one 
insertion of V4, for the case N s = 2. Higher values of N s give similar, but algebraically more 
involved results. In section 2.3 the contribution to the unrenormalized cj) 2-point function from 
those diagrams is denoted fi( 3 ), 

n ab = f ^ \Ao(m 2 a ) + i £ A (m 2 c ) J + ±-5 ab £ (g n a 2 c + g 22 a 2 c ) A (m 2 c ) 

V c=l / c=l 

+ \ {9iiai a a lb + g 22 a 2a a 2b + 3g 12 (a la a 2b + oi 6 a 2o )) ^A (m 2 a ) + A {m 2 b ) + ^Y1 ^o(™c)^ 



, 9 12 1- - - a I 2\ 1 9llll - - \^ ~2 A I 2\ 1 52222 „ * V— <• ~ 2 A I 2\ 

+ —5 ab a 1 -a 2 A (m a ) + — a la a u ^ a i c M m c) + — a 2 a a 2b ^ a 2 c A o( m c) 

c=l " c=l 

9 2 1 2 

+ -^p- ^2 RccabAo{m 2 ) + - {9iii2Slccab + 5l222 S 2cca b) A)(»™c) 
c=l c=l 
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where we defined the rank 4 tensors, 

R = &i (g &i (g a 2 (g a 2 + &i (g a 2 (g &i (g a 2 + &i (g a 2 <8 a 2 (g &i + a 2 (8> &i (g &i (g a 2 

+ a 2 (g &i gi a 2 g> &i + a 2 <g a 2 <g &i <g ai 
Si = a! ax (g ax (g) a 2 + ai <g &i (g a 2 (g a! + hi (g a 2 (g (g a! + a 2 <g hi <g hi g> 

'S'2 = 'S'l|ai^a2 

Also in section 2.3 the contribution to the unrenormalized <p-~x three-point function from diagrams 

(2) 

with one insertion of V4 is denoted (pi,p 2 )- Its explicit form in the case Nf = 2 is the following, 

N f N f 

i<j=l i<j<k<l=l 
1 f N " 

Vab(Pi,P2) = -^ H abB ((pi +p 2 ) 2 ,m 2 a ,ml) - —5 ab J2H cc B ((pi + p 2 ) 2 , m 2 c , m 2 c ), 

c=l 



V i J b 1) ( P i,P2) = -^^a la a lc F 6c i?o((pi + P 2) 2 ,mlm 2 c ) 



6 

c=l 



c,d=l 

v ( a r>( P1 ,P2) = v^>(p 1 ,P2)\ 



— ^2a lb a lc H ac B (( Pl +p 2 ) ,m a ,m c )-—^2a la a lb H cc B (( Pl +p 2 ) ,m c ,m c 

c=l c=l 

— <5 a6 ^ ai c ai d iJ cdJ B ((pi +p 2 ) 2 ,m;?,m^), 



vS 2) (Pi,P2) = -512 < —a la a 2b ^2 H cc B (( Pl + p 2 ) 2 ,m 2 c ,m 2 c ) + — ^ H cb (a la a 2c + a lc a 2a 

{ c=l c=l 

1 Ws 1 

xB ((pi +p 2 ) 2 ,m 2 c ,ml) + —(5 a;) ^ H cd a lc a 2d B ((pi + p 2 ) 2 ,m 2 c ,m 2 d ) \ + (a <-> 6), 
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c,d=l 

V^ 1U) (pi,P2) = -— ai a ai ft ^ a lc ai d H cd Bo(( Pl +p 2 ) 2 ,m 2 c ,m 2 d ), 

c,d=l 

vr 2, (f..p 2 )=v»'"> (P1 ,p 2 )|. . , 

I ai— >a2 

V iI 122) bi'P2) = ^ R abcdH cd Bo(( Pl + p 2 ) 2 ,m 2 c ,m 2 d ), with R defined above, 

c,d=l 

f N " 

V lt- 112) bi'P2) = -— ^ Si ahcd H cd B${(pi +p 2 ) 2 ,ml,m d ), with Si defined above, 

c,d=l 

ai— >a2 
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Gathering the e-poles, we obtain the divergent piece of the H counterterms, 

(8vr 2 e) AH = 1(h+ (H^j + {a^Ha^ ®a 1 ) + 9 -?f^ (fii - a 2 )) 

+^ ( 01 (8) i2a x + i3"ai ® a x + ^Tr (If) a x <g> a x + ^ai-i/ai ] + ^ (&i -► a 2 ) 

+^ 2 - {a\ <8> i2«2 + -H"a 2 (g> a x + ^Tr (if) ai (g) a 2 + ^ai -Ha^j + ^ (oi <-»■ a 2 ) 

.£71112 ,f?1222 

H — — {{ai-Ha 2 + a 2 Hai)ai <g> &i + a r .Hai(ai ® a 2 + a 2 (g) di)) H — (oi a 2 ) 

8 o 

A.l One-loop renormalization of V4(</>) 

In this section we briefly discuss the form of counterterms to V4 in MS scheme and check that those 
counterterms have the same form as V4 at one-loop level. The explicit computation of counterterms 
to V4 in other schemes, such as OS, can be somewhat involved due to their finite parts violating 
Z 2y symmetry. The form of the potential V4 in the classical theory is (2). The potential and its 
counterterms in terms of renormalized fields is given in (10), which can be rewritten as 

V A (z 1 / 2 (t>,{U k a k },v%g + 5g),{iJ%g ij +6g ij )},{v%g ijkl + 5g ij M)}) , (A.2) 

with Uk = Q t Qa k anci , in MS, Q = 1. At one-loop level, (A.2) has the form, 

V4 (z 1/2 0, {U k a k },fi e g, {n € gij}, {n e 9ijkl}) 

+ V 4 (0, {a fc }, ^ e 6g, {n e 8gij )}, {n e 5g ijM )}) + higher orders (A.3) 

The second term has the same form as V4. We want to see that this is also the case for the first 
term, whose explicit expression is, 

N f 



I (0-Z0) 2 + g l3 (fii.C/JzVV) (fij-l/J-Z 1 /^) (4>SZ4>) 

• i<j=i 

+ 1 E (ai-^ 1/2 0) (a j -c/*zv^) (a^z 1 /^) (*vu\z^d>) . (aa) 

The one- loop value for SZ is given in (39) and for SUj = SQ aj m (46). Substituting those 
expressions in (A. 4) we get 

N f 

V A ((f), {a k },n e g, {n e gij}, {^9ijkl}) ~ 4 , 27r 2 6 ^ a ^ &k '^' ^ 

77 e k=i 

1 Nf ( Nf \ 1 Nf 

4!87r 2 e E Si i I E a ™ ( x ™» + ) o-i-^ay^cf) 2 - 4}A/jT2e E fyOfcai -4>a r 4> (a k ■ 4>) 2 

i,j=l \m=l J ' i,j,k=l 

1 N f / N f \ 

~ 4!8^ E ( E a ™ + x ™i + x ™k + x mi) I a* ■ Oj • «fe ■ a r + h.o. , (A.5) 

i,j,k,l=l \m=l j 

where h.o. refers to terms of higher order. (A.5) clearly has the same form as (2). 
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B Loop integrals 

In this appendix we give a list of loop integrals used in the foregoing. More complete calculations 
can be found, e.g., in [3, 7, 8, 9]. Divergent integrals are separated in a dimensional regularization 
pole term and a finite remainder. 71 = //v / 47re~T£ . 



Mm ) = Wf] d V-^ + * £ = + jg^ooCm ), a (m ) = log 



- 1 



B„(p", m ;,m?) = J^L | <f < ( p _ mj - fe) ( i - ?) 2 _ ^ - fe) = + ^f^mf.ml) 

&o(p 2 , m-?, m 2 .) = / dxlog ( (1 — a;)^ 1 + — x(l — s- — ie J 

V M A 4 A 4 / 

W,m 2 ,m 2 ) = J^y d fd d e p _ ml + i£) * + p)2 _ ml + - ) = ^ - ^(/X,™!) 

, 9 2 2\ \ m ? m 2 \P 2 \ 

JO V A 1 A* A 4 / 

which can also be written, 

Br(^,m 2 ,m|)=p^ 1 (/,m 2 ,ml) 
p 2 Bi(p 2 ,mj,ml) = - (A Q (ml) - A (ml) - (p 2 + m\ - m\)Bn{p 2 ,m\,m\)) 

C (pi,P2,m 2 l ,m 2 .,m 2 3 ) = [d d £— , 1 , ■ , 

(2TT) a J (£ 2 -m( + ie){(£ + p 2 ) 2 -m% + ie){(£-pi+p 2 ) 2 -mi + ie) 



C Counterterms to a regular coupling matrix 

The counterterms 5Zjf and <5V to the coupling matrix H, as defined in (12), can be found by 
solving (26) along the lines of appendix C of [1], where the more general case of complex normal 
matrices H is considered. In this appendix we point out that the solution can be simplified if H is 
regular. Our expansion parameters are the (dimensionful) eigenvalues of H. In the basis in which 
H is diagonal, H = diag(/ii, . . . , h^J, we have hj — ► 0, hi/hj = 0(1). We define, 

e = V\hi\ 2 + --- + \h Ns \ 2 = yjTi(H^H) = \\H\\, 
as our perturbation parameter. We want to find SZh from the equation, 

SZ H H +[SV,H]= X, [SZ H ,H] = 0. (C.l) 
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Here, X stands for a known matrix. In the case considered in this paper, X is given by the r.h.s. of 
(26) and it is real symmetric. Since we assume H to be regular, SZh must be a linear combination 
of H k with k = 0, . . . , N s — 1. At one- loop level, SZ H = 0{e 2 ) so we write, 

SZ H = Sz l + S Zl H + 5z 2 H 2 + 0(e 3 ). 

In order to project (C.l), we orthogonalize the basis {1, H , H 2 } of matrices commuting with 
H and being 0{e 2 ). The orthogonal set has the form N = 1, iVi = H - l/N s Tr(H)l, N 2 = 
H 2 —u\H+uq\ j with appropriately determined. When H is normal, Tr H] H n H^ m ) = 

for n,m > 0. Thus, the projections of [5V,i3"] over iVo,i,2 vanish, and we can obtain SZh by 
orthogonal projection from (C.l). 

D Decomposition of a vector function 

Assume we are given a function / : R N — > 1^ of the form /(as) = Ax, where A G R NxN can 
be 03 dependent. We can always find a scalar function A = A(as) > and an orthogonal matrix 
Q = Q{x) G SO(N) such that /(as) = XQx. If the matrix A is close to the identity, A = 1 + <5A, 
to lowest non-trivial order in 8 A we have, 

A = 1 + <5A, 5X = x-SAx, Q = exp(—SQ), 8Q = (6 Ax) A x, 

where x = as/||as||, and u Av = u ® v — v®u. It is clear that 5Xx + SQx = SAx. 

E RG fixed points 

The evolution of Yukawa versors under the RG flow is discussed in section 2.5. In this appendix 
we consider fixed points of the third type, and show that they are saddle points of the flow. The 
beta function for the Gram matrix x is given in (54). From there, we obtain, 

d(3 x .. 

^2_^j_ = 2a 2 ($ u ( x _ XkiXij ) + 6ji(x ki - x kj Xij)) + 2a 2 {8 ik {xij - x u xij) + 5 jk (x H - xijXij)) 
ox k i 

N f 

~^2 a m ( x mi + x mj) ( S ikSjl + S U S jk ) . (E.l) 
m=l 

We consider (3 X as a function of the ( N 2 f ) variables {xij)i<j. Thus, (3 X : —> r( 2 / ) ) and (E.l) 

are the entries of a ( N ^) x (^) matrix. 

At a fixed point of the third type, the Yukawa versors lie on orthogonal directions, with at least 
two of them being collinear (xij = 1 for some i < j), and at least two being orthogonal (x k i = 
for some k < I). We can always relabel {a k } so that . . . , d ni } are collinear, then along an 
orthogonal direction {d ni+ i, . . . , a ni+ri2 } are collinear, and so on. If the {d^} lie along r different 
orthogonal directions in flavor space, n± + . . . + n r = Nf. 

With that convention, the matrix x° has r diagonal blocks of dimension rii x n^, i = 1, . . . , r. In 
each of these diagonal blocks all entries are equal to 1, and outside the diagonal blocks all entries 
are 0. Notice that we have chosen all x°j > 0, as explained in section 2.5. It is easy to check 
that such as satisfies (55). We consider first the case of two orthogonal directions, r = 2. The 
generalization to r > 2 presents no difficulty and is discussed below. 
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In the case r = 2, x° has two diagonal blocks of sizes ni, n 2 with n\ + n 2 = Nf. Assuming 
n i,2 > 1 f° r the moment, from (E.l) we find that d(3ij/dxki has three diagonal blocks, outside of 
which all entries are 0. The first one consists of the entries dPij/dx^i with 1 < i < j < rii and 
1 < k < I < n\. This submatrix of dimensions (™ x ) x (™ x ) is diagonal, with diagonal entries all 

equalto-2E™ 1 =iam<0. 

There is another diagonal block consisting of dPij/dx^i with n\ + 1 < i < j < n\ + n 2 and 
n\ + 1 < k < I < n\ + n 2 , of size (^ 2 ) x (^ 2 ) , which is also diagonal with all diagonal entries equal 
to —2 ^2m=ni+i a m- If "-2 = 1 this block is not present. 

There is, finally, another diagonal block consisting of dfiij/dxki with 1 < i, k < n\ and ri\ + 1 < 
j,l < ni + ri2- This block has size (n\n2) x (niri2) (notice that + (™ 2 ) + n\ri2 = ( N 2 f ))- From 
(E.l), with the explicit form for x° given above, we obtain, 

ni ni+n2 / f>a \ «i+n2 

87f2 E E TT^ =E fl ™- n 1 + l<i<n 1+ n 2 . (E.2) 

k=U=n 1+ l V ° Xfci m=l 

Thus, the sum of the entries of any row of this submatrix is the same for all rows, given by the 
r.h.s. of (E.2) which is independent of (i, j). This implies that we can write this submatrix as, 

(n-i+n-z \ 
"I 1, (E.3) 
m=l / 

where the (niri2) x (71-1^2) matrix A is such that the sum of its columns vanishes. We conclude 
that A has a null eigenvalue and then, from (E.3), the diagonal block we are considering must 
have Ylm=i 2 a m > as eigenvalue. The other eigenvalues of this block can be positive, negative 
or zero, depending on the values of the couplings a^. 

Since we are considering a fixed point x° of third type, at least one of n±, n 2 must be > 1, and 
we can always assume that n\ > 1. From the foregoing, we see that independently of the values 
of the coupling constants a^, and of whether n 2 = 1 or n 2 > 1, the linearized beta function has at 
least one negative eigenvalue —2 Ylm=i a mi an d at least one positive eigenvalue Y^=\ 2 a rn- Thus, 
these fixed points are saddle points. 

The generalization of this conclusion to the case r > 2 is immediate. We have diagonal blocks of 
size ("*) x C£) with eigenvalues —2 Ylm^ni+^+ni^+i a m > 1, and blocks of size (ninj)x(ninj), 

i / j, with one eigenvalue equal to Z)m=ni+"!+Jii_i+i a m "I" E m =ni+.. i +nj_i+i a m- (-^Y convention 
we set no = 1, if it occurs in these expressions.) In this general case, too, the linearized beta 
function has negative and positive eigenvalues at a fixed point of third type. 



F The su(N) algebra 

In this appendix we gather our notations and conventions for the su(N) algebra, and several 
relations among its constants that are used in section 3. We denote k = k 2 — 1 for brevity. 
Gell-Mann matrices for su(N) are denoted A^, 1 < a < N, with the supraindex indicating the 
dimension usually omitted. The diagonal matrices are 



A ^ = \ Hf TT diag(l,...,l,-(m-l),0,...,0), m = 2,...,N. 

m 
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The definition we give for A^ is recursive. A^ 3 are Pauli matrices. Given the set {Aa Vj }^ =1 , we 



,(2) 



obtain with a < N by embedding Aa v; in the upper-left corner of a null (N + 1) x (N + 1) 

matrix. The matrix A^t- 1 ^ is diagonal, as defined above. Finally, the remaining 2N matrices 

\ a N+1) with N <a< WTT are, 



AN) 



X (iV+l) _ 
N 2 ~ 



/o 




\i 



1\ 



0/ 



We have the usual relations, 



y N 2 +l 



/o 




\i 



-i\ 



0/ 



(N+l) 
N+l-2 



/o . 


. 0\ 




/o . 


. 


o\ 


. 


. 


' 'V+i-i 


. 


. 





. 


. 1 




. 


. 


—i 


\0 . 


. 1 0) 




\0 . 


. i 


0/ 



N 



A a A a 



N 2 -l 



N 



N 



o=l 



N 



1> ^ A a AfeA a = — — Aft, y^(Aq) ij(X a )ki = 2$ii$kj 



a=l a=l 

Tr (A a Aft) = 25 ab , Tr (A a A 6 A c ) = 2{d abc + if abc ). 



(F.l) 



The totally antisymmetric structure constants f abc and the totally symmetric constants d abc are 
defined relative to this basis by, 



N 



(d a bc + ifabc)X c + JjSabl 



c=l 



(F.2) 



Particularly important for the evaluation of radiative corrections in section 3 are the constants f abc 
and d abc with one index equal to N. Those are given by, 



( I N 

Y 2(AT-1 



Nab 



UN— 1 < a < N, b = a + 1, and A a is real, 



N 



2{N 



— UN— 1 < a < N, b = a — 1, and A a is imaginary, 



otherwise, 



(F.3) 



and d Nab = d NaJab with, 



Naa L 



dr 



Naa 



w 



2 

N(N-l) 
N-2 
~y/2N(N-l) 

-(N-2 



if 1< a < N - 1, 



N(N-l) 



if TV - 1< a < TV , 
if a = JV. 



(F.4) 



With these equations we can compute the algebraic expressions appearing in loop diagrams in 
section 3. We quote here the results for the divergent parts. In fermion tadpoles we have, with 
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given in (68), 

J2( x Nh m i = 4 JV jv(jv-i) J ( ^ (o))3 ' (R5) 



, N 2 - 3N + 3 



and in fermion bubble diagrams, 

N 



= (^ (0 )) 2 5 a6 2 (1 + (^ aa ) 2 + Mjj^djrjrjj - (fjj a{a+1) ) 2 - (fm(a-i)) 2 + f <w) (F-6) 



(ff«(0)) <U 



1 



iV(iV - 1) 



x < 



12 if 1 < a < N - 1 

4(iV 2 -3Af + 3) ifiV^T<a<iV. 



^12(^2 -3iV + 3) if a = AT 
In boson tadpole diagrams we find expressions of the form, 

f d - (d + 2d - )m 2 - 2 (iV-2)(iV-4) (iV - 2)2 

/ , Ncc\ NNN ^ za Ncc) m c ~ z — 1 N(N — 1) ^ ' 

with m 2 , from (67). In boson bubble diagrams we find expressions of the form, 
17 N 2 -4 

Yl d Naa d Nbb d acddbcd = { d N aaf 5 ab , (F.8) 



iV 

c,d=l 



Y d Naa d Ndd d acddbcd = ( d Nacf 6 CLb, (F.9) 



2iV 

c,c(=l 



iV — 4 

E ^(^W + 2d Ncc) = 2 ^^ S Nb> ( F -!0) 



AT 

c=l 
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E ( d Ncc) 2d acddbcd = Sab ( 1 ~ "T^ + ^ a ]y + TF ( (/jVa(a+l)) + (/jVa(a-l)) ) 



2 /iV _ 6_ , 



_/V JVaa \ 4 iV 

N 



^2 3 



'Xnu 1 N d Naa d NNN 



(F.12) 



E d Ncc d Ndd d acddbcd = Sab (o ~ ^ + \ T ~ Jj) ( d N, 

■4=1 V V 7 

I ~ iV") ((^(«+l)) + (-^Va(a-l)) J ^ \2~ WJ" Na J ' 

Y Y d abl d nnl = - Nd TMN d Nab> ( F - 13 ) 
n=l i=l 

/V_1 ^ AT 3A 7 " — 2 JV_1 N — 2 

Y E = -Y d 7rWN d abN ~ N (N_i) Sai > Y 6an + JfZT[ SabS <ar + ( N ~ 2 ) 6 ^- (F.14) 

n=l J=l ^ ' n=l 



JV-1 iV 
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An expression related to the last one is, 
37=1 W N-2 I 2N 1 37=1 7V-2 

/an*/&n* = ~^—]J ]y37<W + [<U <U - -y^b^ + ( N ~ 2 )<U- (F.15) 

n=l i=l ' n=l 

All of these relations are easily checked numerically for low values of N. 
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